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Abstract 

The states of asymptotic relaxation of 2-dimensional fluids and plasma present 
a high degree of regularity and obey to the sinh-Poisson equation. We find that 
embedding the classical fluid description into a field-theoretical framework, the same 
equation appears as a manifestation of the self-duality. 

The states generated by externally driving (stirring) an ideal fluid can have very 
irregular form. It is however known from experiments and numerical simulations that after 
suppressing the drive the system evolves to states with high degree of order, essentially 
consisting of few large vortices, with very regular geometry. These states are attained 
after long time evolution and are not due to the residual dissipation. The process consists 
of vortex merging, which is an essentially topological event where the weak dissipation 
only allows the reconnection of the field lines but does not produce significant energy 
loss from the fluid motion. Inferring from results of numerical simulations, Montgomery 
et al. JIJ, @ have proved that the scalar stream function ip describing the motion in 
two-dimensional space obeys in the far asymtotic regime (where the regular structures 
are dominant) the sm/i-Poisson equation 

Aifj + 7 sinh (f3vb) = (1) 

where 7 and (3 are positive constants. The relations of ip to the velocity and vorticity 
are v = xe z ,u; = Vxv = —V 2 ipe z where e z is the unitary vector perpendicular 
to the plane. With these variables, the Euler equations for the two dimensional ideal 
incompressible fluid are 

_ + ( v . V )u, = (2) 
V-v = 



We will try to develop a field-theoretical model of the stationary asymptotic relaxed 
states, i.e. we look for a model that could provide a derivation of Eq. ([[]). This will be done 
progressively, examining models ellaborated for closely related problems and collecting the 
relevant suggestions that could allow us to write a Lagrangian density. 

In the study of the two-dimensional Euler fluids, and in particular in explaining the 
origin of Eq.(|l|), an important model consists of a system of iV discrete vorticity filaments 
perpendicular on plane, having circular transversal section of radius a and carrying the 
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vorticity Ui, i = 1, N . A very comprehensive account of this system is given by Kraichnan 
and Montgomery |§ where the correspondence between the continuous and discrete rep- 
resentations of vorticity is discussed in detail. The motion in plane of the k-th filament 
of coordinates r k = {r\, r%) = (x^, i/k) is given by 

-dt = £i3 7fl E ^G(r fc -r n ) , i,j = 1,2, k = l,N (3) 

where the summation is over all the other filaments' positions r n , n ^ k, and e*- 7 is 
the antisymmetric tensor in two dimensions. As shown in Ref . || G (rjt — r n ) can be 
approximated for a small compared to the space extension of the fluid, L, o < I, as the 
Green function of the Laplacian 

GforOw-^ln/t^T) (4) 



2i \ I 

The statistical properties of the system of discrete vortices have been examined using 
the Liouville theorem and the conservation of energy and momentum . The model consists 
of an equal number of positive and negative vortices with equal absolute magnitudes |u^| = 
in contact with a thermal bath of temperature T. The possibility of the generation of 
two supervortices of opposite signs arises when T is negative and sufficiently large. When 
the most probable state is attained for a stationary configuration the stream function if) 
is shown to verify the sinh-Poisson equation ([!]). These statistical considerations remain 
the reference explanation for the appearence of this equation in this context. 

Highly ordered states have been obtained by numerical simulations for other systems 
in two spatial dimensions: systems of filaments of electric current in 2-dimensional MHD, 
guiding centre particles. The asymptotic states are found to be described by scalar func- 
tions obeying Liouville equation or sinh-Poisson equation. The theoretical approach which 
permitted to explain this result was based on the extremum of the entropy over a statis- 
tical ensemble of states of a set of discrete objects under the constraints which express 
the conservation of the invariants of motion and with a negative temperature. 

In the equations of motion (^) the right hand side contains the curl of the Laplacian 
Green's function (f|) (we take L — 1) 

-e%GM=e*d j ±]nr = ±e«£ (5) 

V 2 — lnr = 5 2 (r) 

2tt v ; 

The term in right hand side of (J3|) can be seen as a vector potential a (r, t). The "magnetic" 
field V x a is a sum of Dirac 5 functions at the locations of the vortices. If we take equal 
strength uo for all vortices this potential appears as the "statistical potential" and has a 
topological interpretation ||. From fl5|) it can be rewritten 



1 • • r J 1 y 1 

— £ v — = ftarctan— = dS (6) 

2tt r 2 2tt x 2n 
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where r = (x,y) = (r cos 9, r sin 6). The "magnetic" flux through a surface limitted by a 
large circle is proportional with the number of vortices. The topological nature of this 
potential suggests that it can be naturally derived in a topological framework, i.e. from 
a Lagrangian density of the Chern-Simons type, 

C = ~e^A w F a0 (7) 

where F a p = d a Ap — dpA a . In the analysis for the two-dimensional Euler equation (as 
in the investigation of similar systems), the vorticity-type functions uj are expressed as 
the Lapacian operator applied on a scalar function and this naturally invokes the Green 
function of the Laplacian. Then this approach suggests that the intrinsically determined 
motion of the fluid can be projected onto two distinct parts of a new model: objects 
(vortices) and interaction between them (the potential obtained from the Green function). 
The fact that the interaction potential can be derived from a Chern-Simons topological 
action suggests that, in order to embedd the original Euler fluid system into a larger 
field theoretical context, we need (a) a "matter" part in the Lagrangian, which should 
provide the free dynamics of the vortices; (b) the Chern-Simons term, to describe the 
free dynamics of the field; (c) the interaction term of the (Chern-Simons-) field and the 
matter. The main requirement to this field-theoretical extension of the original Euler fluid 
model is to reproduce the discrete vortices and their equation of motion. 

Jackiw and Pi ||, || have examined a model of N interacting particles (of charges e s ) 
moving in plane described by the Lagrangian 

N i i r r 

where m s v s = p s - e s A (r s \r 1 , r 2 , r N ), A (r s |r a , r 2 , r N ) = (a* (r 1 , r 2 , riv) i=lj2 and 
a* (ri, r 2 , rjv) = ^£ lj J2qjt s e ? jT^ZT^- ^he matter current is defined as j M = (p,j) = 

J^fLi e sV s 5 (r — r a ), v% = (1, v a ) with the metric (1, —1, —1). By varying the action we 
obtain the following equations 

^F af3 = e^d a A p = f 

B = -p (9) 

E i = e ij j j (10) 

Here B = e^djAi. In these equations the potential appears as being generated by the 
matter, i.e. by the "current" of particles. We note that the Eq.(|9]) connects the matter 
density with the curl of the potential. This is important since in order to derive it in the 
context of the quantum version of their model, Jackiw and Pi have shown that one needs 
to include a nonlinear self-interaction of the wave function representing the matter field. 
The self-interaction of the scalar matter field is of the type <£> 4 . The canonical momentum 
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expression becomes the covariant derivative in the field theory version. Even at this point 
where we only have some hints, it is suggestive to look for possible identifications between 
the fluid variables and the field theoretical variables of the model of Jackiw and Pi. The 
fluid variables are: the scalar potential i/j, the velocity v l = e^dji/j and the vorticity 
uj = e^diVj. In order to check the possibility that the fluid model can be embedded into 
the field theory just described we identify 

u = p = ^ (11) 



v % = A' (12) 

Consider the relation between the vector potential A and the density p, obtained from 
Eq.© 

tfdiA' = -e lk d kP 

or, writting symbolically the Green function of the Laplace operator for argument r — r' 
as the inverse of the operator at the left, we have 



Ai (r, t) = e> k d k J d 2 r' (-8%)^, p (r', t) 



(13) 



Since uj = — V 2 ^ (where V 2 = d l di) we have formally if) = — (V 2 ) 1 uj and use this 
formula to calculate the right side term of (13) taking into account the identification flllD 



^4 /«"(-**)><*) = .»<W 



V 1 



which confirms (12). It results that if we assume the identification (11) the equation (|9]) 
obtained form the Lagrangian is precisely the definition of the vorticity vector. 

One important hint from the model of Jackiw and Pi is the idea to represent a classical 
quantity as the modulus of a fictitious complex scalar field, as in Eq.(|ll]) and derive 
dynamical equations from a Lagrangian density expressed in terms of this field. The 
substitution of the dynamics expressed in terms of usual mechanical quantities by the 
richer dynamics of the amplitude and phase of the complex scalar field is an example of 
embedding of one theory into a larger framework and relies on the example of quantum 
mechanics. 

In Ref.(@) the sinh-Poisson equation is derived in a field theoretical model where 
instead of the topological coupling two complex scalar fields are considered. In that 
model the dynamics of the two scalar fields is independent except that the self-interaction 
depends on both. Combining this suggestion with the one exposed in the previous para- 
graph, we will take the density p (which we identify with the local value of the vorticity 
uj as in Eq. (|Tl|) ) as a two-component complex field ^ 

p ~ [Tjrt ) ^rj 
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where = J 1 ) . Keeping the same structure of the Lagrangian density as in the 

V 02 J 

previous example but extending to a non-Abelian SU (2) gauge field we have ||, |J 
C = -e^tr (d,A u A p + -A^A/J + itr (^D V) - Ur ((A^) f A*) + \tr([&, 

(14) 

where the potantial A takes values in the algebra of the group SU (2). This form in- 
corporates and adapt all the suggestions from the previous models: (a) the first term is 
the general non-Abelian expression of the Chern-Simons term; (b) it uses the covariant 
derivatives (/x = 0, 1, 2, i = 1,2) for the minimal coupling 

D^ = d^ + [A^] 

(c) Finally, it includes a scalar self-interaction of ip 4 type. Analogous to the case treated 
by Jackiw and Pi for the U (1) gauge field @, in Ref. || it is shown that the Hamiltonian 
density corresponding to the Lagrangian density flU]) is 

H = \tr ((A^) f (A*)) - \tr ^] 2 ) (15) 

since the Chern-Simons term does not contribute to the energy density (being first order 
in the time derivatives). The equations of motion are 

iD^ = - l -VH- l -[[^!^]^] (16) 



Ffj,u — c^iivpJ 9 (17) 
Using the notation D± = A ± 1D2 the first term in ([T5) can be written 

tr (( A^) f (A*)) = tr ((.D-^) 1 + [ [*> . *] ) 

The last term comes from Eq.fllTD and from the definition J° = ftf, . Then the energy 
density is 

ft = itr ((£>-^) f (!>-*)) > (18) 
and the Bogomol'nyi inequality is saturated at self-duality 

D^ = (19) 

d + A_ - d-A + + [A+, A-] = [*, (20) 
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The first equation results from the minimum in Eq.(^) and the second is actually Eq.fll 
with the definitions 



J° 
J* 



[*t,*] 



The static solutions of the self-duality equations C pLSj , p0|) are derived in ||, using the 
algebraic ansatz: 



a=l 



M 



a=l 



where H a are the Cartan subalgebra generators for the gauge Lie algebra, E a are the 
simple-root step operators and E_m is the step operator corresponding to minus the 



maximal root [10] . The rank of the algebra is noted r, and r = 1 for SU (2). Then 



-M 



(21) 



a=l 



The equations (19, 20) lead to the affine Toda equations 



r+l 



(22) 



6=1 



for a = 1, r, plus the index for M, i.e. a = 1, 2. C a b is the extended Cartan matrix 

2 a («) . a W 



o6 



(6) 



|2 ' 



a,b= 1,2 



where a;'-' 1 -' are the simple root vectors of the algebra sit (2), and in addition the minus 

r | 2 

(23) 

Alnp 2 + 2(- Pl + p 2 ) = 
and this gives the relation 



maximal root. The equations ( |22| ) can be written in detail for p 1 = l^ 1 ) 2 , P2 = \i> AI \ 2 

Aln Pl + 2( Pl -p 2 ) = 
Alnp 2 + 2(- Pl + p 2 ) = 



A In ( Pl p 2 ) = 



or 



p 2 = const p 1 1 



(24) 
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since the exponential of a linear term is excluded by the conditions on a circle at infinity. 
Using Eq.(^l|) we have to identify 

w = Pi ~ P2 (25) 
whose equation is obtained from equations in (|23| ) using ( p4|) const= 1 (see bellow) 

A In pi + 2( Pl -pf 1 ) = (26) 



The substitution ^ = lnp x leads both equations (^5|) and ( p6[) to the form 



A^ + psinh (|/?|^) = (27) 



Choosing a constant different of unity in Eq. (p24|) imposes a linear substitution of the 
stream function : ifi — > if)' = 7 + and modifies the factor multiplying the sinh function 
in Eq.flTTp. 

We have constructed the Lagrangian density (fill) with a standard structure of matter- 
gauge field interaction and incorporating suggestions from models relevant for our objec- 
tive. The aim was to reproduce the model of fluid vortex filaments supposed to be a 
faithfull representation of the original Euler fluid dynamics. The fundamental step in the 
derivation of the sinh-Poisson equation is the assumption that the energy density (|18|) is 
minimum (which is equivalent to minimizing the action for these stationary solution) i. e. 
in the condition of self- duality. 

We conclude that the Euler fluid equations can be embedded into a field theoretical 
framework via the discrete model of vortex filaments. This framework is able to provide 
the equation obeyed by the scalar stream function at stationarity and this is precisely the 
equation inferred from numerical experiments. The particularity of the field theoretical 
framework is that it shows that the asymptotic states consisting of regular vortices corre- 
spond to self-dual states in field theory. This might have deep significance and deserves 
further investigation. 

One must also note that the statistical-theoretical model able to explain the sinh- 
Poisson equation [pl| , [|12T] , |fL3| , Ulij appears to have strong and deep connections with 



the self-duality. It should be carefully considered as a possible approach in the study of 
models where self-duality is known to exist, like O (n) or Self-Dual Yang-Mills theory. 
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